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Abstract 

We consider the nonlinear Dirac equations in one dimension and review various 
results on global existence of solutions in H . Depending on the character of the non- 
linear terms, existence of the large-norm solutions can be extended for all times. Global 
existence of the small-norm solutions is proved for the most general nonlinear Dirac 
equations with cubic and higher-order nonlinear terms. Integrability of the massive 
Thirring model is used to find conditions that no solitons occur in the Cauchy problem 
with small initial data in a subspace of L 2 . 

1 Introduction 

The goal of this article is to survey a number of recent results on global well-posedness of 
the nonlinear Dirac equations in the space of one dimension. The nonlinear Dirac equations 
are known for long time in quantum mechanics and relativity theory [22]. Recently 
these equations were used to model other physical systems such as photonic crystals and 
Bose-Einstein condensates in optical lattices [TJ [T71 [2T] . 

The nonlinear Dirac equations are similar to the nonlinear Klein-Gordon equation on 
one hand and to the nonlinear Schrodinger equation on the other hand. In the former case, 
the reduction is possible for a special form of nonlinear terms as the linear dispersion rela- 
tion between the two models are identical. In the latter case, the reduction appears in the 
asymptotic limit using an envelope wave approximation near one branch of the wave spec- 
trum. Since analysis of global well-posedness in the nonlinear Klein-Gordon and nonlinear 
Schrodinger equations has been booming in the last ten years, it is not surprising that the 
interest of harmonic and PDE analysts turns recently to the nonlinear Dirac equations. 

The organization of this article is as follows. Section 2 sets up the nonlinear Dirac 
equations and reviews a number of physically relevant models. Global well-posedness in H 1 
is studied in Section 3 using apriori estimates as in Goodman et al. [6]. We show that the 
nonlinear Dirac equations with a special structure of nonlinear terms are globally well-posed 
even for large powers of the nonlinear terms. This result is different from the behavior of 
the nonlinear Schrodinger equation. 
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Section 4 deals with global solutions for small initial data in H using analysis in 
Strichartz spaces as in Pelinovsky & Stefanov [18] . We prove the global existence of small- 
norm solutions that scatter to zero for the nonlinear Dirac equations with quintic and 
higher-order nonlinear terms. 

The nonlinear Dirac equations with cubic terms are considered in Section 5. Scattering 
to zero for small-norm solutions is proved using the recent results of Hayashi & Naumkin 
[SJ IS]. We also discuss the special role of the integrable version of the nonlinear Dirac 
equations known as the massive Thirring model [22], for which global well-posedness in L 2 
was proved recently by Candy [2J. Using the formalism of the inverse scattering transform, 
we show that small initial data in a subspace of L 2 are associated to purely continuous 
spectrum of the Lax operator and admits no solitons in the long-time asymptotics. 

Acknowledgments: Part of this work was done during the visits of the author to 
University of British Columbia and the Kyoto University in 2010. The author would like 
to thank K. Nakanishi and T.P. Tsai for useful discussions regarding this article. 



2 Model 

Let us consider the nonlinear Dirac equations 



i(u t + u x ) + v = d u W(u, v), 
i(vt - v x ) + u = dj}W(u,v), 



(2.1) 



where (x,t) G R 2 , (u,v) G C 2 , and W(u,v) : C 2 
the following three conditions: 



is a nonlinear function which satisfies 



• symmetry W(u,v) = W(v,u); 

• gauge invariance W(e l9 u,e ie v) = W(u,v) for any 9 G M; 

• polynomial in (u,v) and (u,v). 

The nonlinear Dirac equations can be rewritten in the abstract evolutionary form, 



idtu = Hu + f (u) 



H 



-id~ 



-1 



■1 id x 



u 



u 

V 



f(u) 



d a W(u,v) 
djjW(u, v) 



(2.2) 



A homogeneous quartic polynomial W(u,v) satisfying the three properties above is 
characterized by Chugunova & Pelinovsky [3], 



W = ai(|n| 4 + \v | 4 ) + ct2 |w| 2 1 v | 2 + as(uv + uv) 2 + a^v^ 2 + \v\ 2 ){uv + uv), 

where (a±, ct-ii a 3> a i) G ^ are numerical coefficients. 

The standard example of the nonlinear Dirac equations with 



(2.3) 



W = a(\u\ 2 + \v\ 2 ) 2 + 2a\u\ 2 \v\ 2 , aG 



(2.4) 
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occurs in the context of periodic dielectric materials under the Bragg resonance [21] . The ac- 
count of zero-mean periodic modulation of the nonlinear refractive index gives the nonlinear 
Dirac equations with 



W = a(uv + uv)(\u\ 2 + \v\ 2 ) + f3((uv + uvf - 2|n| 2 |u| 2 ), 



(2.5) 



where a E R and j3 E M are proportional to different Fourier coefficients of the nonlinear 
refractive index [lj. 

Two models are relevant for general relativity: the massive Thirring model when W = 
4|ii| 2 |-u| 2 [22J and the massive Gross-Neveu model when W = 2(uv + uv) 2 [TJ. Under the 
following change of variables, 
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iP = Tm, g = Tf , T 
the nonlinear Dirac equation (|2.2p can be written in the equivalent form 

id t tp = Mip> + g, M = 



1 d x 
-d x -1 



(2.6) 



(2.7) 



If if} = ("0, (f>), system (j2.7|) can be written as follows: the massive Thirring model is 

iip t - ip - <p x = (ip 2 + (j) 2 )^, 
i(f) t + <f> + ip x = (ip 2 + <p 2 )^, 

and the massive Gross-Neveu model is 

iipt-rp- 4> x = (ip 2 - <j) 2 )ip, 
i(p t + (p + ip x = {(p 2 - ip 2 )cp, 



(2- 



(2.9) 



Thus, we see that all coefficients of the general quartic potential (12. 3p have physically 
relevant applications. Note also that W may also have sixth-order and higher-order terms, 
as in the context of the Feshbach resonance for Bose-Einstein condensates [19], where 



W = a(\u\ 2 + M 2 )|u| 2 M 2 , 



Q E 



Local existence of solutions of the nonlinear Dirac system in Sobolev space fP(R) can 
be proved with standard methods using the Duhamel formulation and the fixed-point ar 
guments [HE]. If uo E H s ( 

nonlinear Dirac equations (|2.2p admits a unique solution 



for a fixed s > |, then there exists a T > such that the 



s-l 



u(t) E C{[0,T],H s (R))nC l {[0,T],H 

where u(i) depends continuously on the initial data u(0) = uo- 

In what follows, we review results on global well-posedness of general nonlinear Dirac 
equations (|2.ip in some subspaces of if^R) or L 2 (R). 
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3 Global well-posedness in H 1 

There exist three conserved quantities of the nonlinear Dirac equations (|2.ip corresponding 
to Hamiltonian H, momentum P, and charge Q, 

H = o / {u x u — uu x — v x v + vv x ) dx + / (vu + uv — W(u, v)) dx, (3.1) 
2 Jr Jr 

P = — (uu x — u x u + vv x — v x v) dx, (3.2) 
2 Jr 

and 

Q= f (\u\ 2 + \v\ 2 ) dx. (3.3) 
Jr 

These conserved quantities are well defined for a local solution in iiP(R) for s > \ thanks 
to the Banach algebra of H S (M.), s > i with respect to multiplication. 

Unlike the nonlinear Schrodinger equation, the Hamiltonian H is not useful for analysis 
of global well-posedness because the quadratic part of H is sign-indefinite. Nevertheless, 
for a special nonlinear function W(u,v), local solutions can be extended to global solutions 
for all t £ 1. The following theorem generalizes the result of Delgado [4] for the massive 
Thirring model (|2.8|) and the result of Goodman et al. [6] for the nonlinear Dirac equations 
with W in (El). 



Theorem 3.1 Assume that W is a polynomial in variables \u\ 2 and \v\ 2 . Let u(t) £ 
C([0.T], H 1 (R)) be a local solution of the nonlinear Dirac equations 112. 1\) for some T > 0. 
Then, u(t) 6 C(R, H l ( 



Proof. To extend the local solution u(t) G C([0, T], il^R)) for all f e R, it is sufficient 
to prove that the H 1 -norm of the solution u(t) satisfies the estimate 

sup ||u(i)|| H i < C(T), (3.4) 
te[o,T] 

where the constant C(T) is finite for T < oo but may grow as T — )• oo. By the conservation 
of Q in (13. 3p . we have 

||u(t)|| i2 = ||u |U2, teR. (3.5) 

To consider ||u(i)[|^2 P +a for a fixed p > 0, we multiply the first equation of system (|2.ip 
by |ii| 2p u and the second equation by \v \ 2p v, add the two equations, and take the imaginary 
part. If W depends only on \u\ 2 and \v\ 2 , the nonlinear function is cancelled out and we 
obtain 

_J_$ (| n |2p+2 + ^12^+2) + J_g x n u f P+ 2 _ | v |2p+ 2) = ./ - _ ^ )( | n |2 P _ ^2py 

Integrating this balance equation on x G R for a local solution in C([0, T], fl" (R)) and using 
inequality \u\\v\ < ^(|«| 2 + |v| 2 ), we obtain apriori estimate 

|l|u(t)||?+ + 2 2 < 4(p+ l)||u(t)||f + + 2 2 . (3.6) 
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By Gronwall's inequality, we have 

\\u(t)\\ L 2 P +2 < e 2|t| ||u || L 2 P+ 2, t G [0,T]. (3.7) 

Since the estimate holds for any p > 0, it holds for p — > oo and gives apriori estimate on 
the L°°-norm of the local solution u(t). The bound on the L°°-norm is needed to control 
the growth rate of the L 2 -norm of u x (t) as t — > oo. 

Taking ^-derivatives and performing a similar computation, we obtain the balance equa- 
tion 

dt {\u x \ 2 + \v x \ 2 ) + d x (\u x \ 2 - \v x \ 2 ) = i{u x d x du + v x d x d^ - u x d x d u - v x d x d v ) W(u,v). 

Let TV" > 1 be the degree of the polynomial W in variables \u\ 2 and \v\ 2 . Integrating the 
previous equation over i£R and using bound (|3.7p . we obtain the estimate 



|||u,(t)||! 2 <C w e^- 1 )l t l||u,(t)||I 2 , 

where the constant Cyy > depends on the coefficients of the polynomial W. By Gronwall's 
inequality again, we obtain 

% f„4fJV-lHtl il „ 

\\Mt)\\h < eTO( ] \\(PoMh, t e [o^i- I 3 - 8 ) 

The exponential factor remains bounded for any finite time T > 0. Therefore, C(T) < oo 
if T < oo and bound fj3.4f> gives global well-posedness of the nonlinear Dirac equations in 
the if 1 -norm. □ 

Remark 3.2 The result of Theorem \3.1\ is very different from the behavior of the nonlinear 
Schrddinger equations, where global solutions may not exist for nonlinear terms generated 
by a polynomial W in variables \u\ 2 and \ v\ 2 of degree N > 3. 

If W also depends on (uv + uv), apriori estimates of the L p -norm include nonlinear 
terms which may lead to the finite-time blow-up of solutions in L°° and H l norms. That 
is, there may exist T max < oo such that 

lim ||u||#i = oo. (3-9) 

The next sections address the question if the finite-time blow-up (j3.9[) can be excluded at 
least for small- norm solutions in the general nonlinear Dirac equations (12. ip . 

4 Global well-posedness in Strichartz spaces 

Strichartz spaces L v t L x and L q x L v t are defined for 1 < p, q < oo by the norms 

T \VP / r \ 1/q 



(J 11/(^)11^^) , \\f\\Lli*-= U\\f{x,-)\\l»dx^ \ (4.1) 
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where T > is an arbitrary time including T = oo. Strichartz spaces have become popular 
in the context of global scattering and asymptotic stability of solitary waves. Nakanishi 
[TO] applied these spaces to the proof of global well-posedness for small initial data in the 
nonlinear Klein-Gordon and Schrodinger equations. Here we consider the nonlinear Dirac 
equations and modify arguments from a more general work of Pelinovsky & Stefanov |18j 
on the asymptotic stability of small solitary waves. 

Let H. be the one-dimensional Dirac operator in (12. 2p and R-^(X) = (H — A/) -1 be the 
resolvent operator, defined as a bounded operator from L 2 (M) to L 2 (M) for any A ^ cr(H), 
where 

<r(H) = (-oo,-l] U [l,oo). 

Using Fourier transform, 

/(*) = 7T- I f(k)e ikx dk, />) = [ f(x)e- ikx dx, 



2vr .. 

the resolvent operator i?%(A) can be expressed in the explicit form 



(RuWf)(x) 



1 f e 



ikx 



2tt L A 2 - 1 - k 2 



-(k + X) 1 
1 (k — A) 



f(x)dx, X $ a(H). (4.2) 



Let k G C be a solution of the algebraic equation /t 2 + A 2 = 1 for A ^ cr(%) such 
that Re(ft) > 0. After computations of the Fourier integrals, the resolvent operator can be 
expressed in the Green's function form, 



(RnWfKx) = 1- ! 

2k Jr 



X + i/«sign(x — y) —1 

—1 —A — iKsign(x — y) 



e 



-^-y\f{y)dy. (4.3) 



The Green's function representation (|4,3p shows that the resolvent operator can be 
extended to the continuous spectrum a(T-L) as a bounded operator from L 1 (M) to L°°(R) 
for any A € cr('H)\{l, —1} excluding the end points ±1 of the continuous spectrum, 

H£(\) :=HmR H (X±ie), X G (-00, -1) U (1, 00). 

Using the same representation (14. 3h . we can see that for any Ao > 1, there is C(Ao) > 
such that 

sup \\R± (A) || < C(A ). (4.4) 

|A|>A 

Dispersive decay estimates for the semi-group associated to the one-dimensional Schrodinger 
operator [23j extends to those e~ ttH . Indeed, both semi-groups have the same dispersion for 
small Fourier wave numbers k. For large Fourier wave numbers, the semi-group e~ lt ^ be- 
haves similar to that of the wave equation and it is controlled using the uniform asymptotic 
behavior (|4,4p . As a result, there is C > such that 



-an 



^Ct-^Wfh,. (4.5) 
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Interpolating with the conservation law ||e **^/|| £ 2 = II/IIl 2 > we obtain for any 2 < p < oo, 



1 1 



f\\rS < Ct^-^WfU,, - + -7 = 1- (4.6) 



p p 



Pointwise dispersive decay estimates (14.6P allow us to introduce Strichartz admissible 
pairs and Strichartz estimates for the nonlinear Dirac equations. 

Definition 4.1 We say that a pair (q,r) is Strichartz admissible for the nonlinear Dirac 
equations if 

2 11 

q > 2, r > 2 and - + -<-. 

q r 2 

In particular, (q,r) = (4, oo) and (q,r) = (oo,2) are end-point Strichartz pairs. 

Lemma 4.2 Let (q,r) be a Strichartz admissible pair. There are constants C > such that 

\\e- itn f\\ LtL ~<C\\f\\ Hl , (4.7) 
\\e- im f\\ L rHl<C\\f\\m, (4.8) 



[ l e~ l(t ~ T ^ H F(T,-)dT 
Jo 



<C||F|| L i^. (4.9) 



The following theorem simplifies the main result of Pelinovsky & Stefanov |18| to the 
global small-norm solutions of the nonlinear Dirac equations with quintic and higher-order 
nonlinear terms. 



Theorem 4.3 Consider the nonlinear Dirac equations 112. 2\) with 

f(au) = a 2n+1 f(u), a G E, 

for a fixed integer n > 2. Assume that u(0) G ff 1 (lR) and ||u(0)||#i is sufficiently small. 
The nonlinear Dirac equations \2. 2\) admits a global solution 



u(t) G C(M+ , i? 1 (M) ) n L 4 (M+ , L c 

Proof. By Duhamel's principle, we can rewrite the Cauchy problem for the nonlinear 
Dirac equations (12. 2p in the integral form, 



u(t) = e" im u(0) + f e-^Wf (u(s))ds. (4.10) 
J o 

By Lemma 14.21 solutions of the integral equation (|4.10p satisfy the bound, 

H u llL4 L ^nL-//i < C\\u \\ H i + C||f(u)|| £ i H i, (4.11) 
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for some C > 0. We set up the problem of solving the integral equation (|4.10p as an iteration 
scheme, where we look for a fixed point in a small ball in normed space LfL^ n LfH].. 
Because f(u) is a homogeneous polynomial of degree 2n + 1, we obtain 



l|f(u)llii/ri < C||(|u| + \d x n\)H n \\ LlLl < C\\u\\ LTHl \\urX Loo . 
By Sobolev embedding and the log convexity of the L 2n (]R) norms for any n > 2, we have 
II II ^ II n 2 / n ii i|i~ 2 / n n\\ II 

As a result, we obtain 

||f(u)||^ < C\\*f^ LTH ^ 
and the fixed point argument is closed for small u(0) € □ 

Remark 4.4 Because ||u(i) a continuous function of t € M+ and ||u(t) ||z,°° £ 
L 4 (M + ), we /icwe 

lim ||u(t) ||x,oo = 0. 

t— too 

Moreover, ||u(i)||.z,oa = 0(i _1 / 4 ~^) as t oo /or some v > 0. 



5 Scattering around zero for cubic Dirac equations 

Here we consider the nonlinear Dirac equations (|2.ip with the cubic nonlinear terms, which 
are generated by the quartic function W in (|2.3p . If W is a function of \u\ 2 and |-u| 2 , results 
of Theorem 13.11 show that the local solutions in H 1 are globally well-posed. On the other 
hand, decay to zero of small initial data is excluded from results of Theorem 14.31 because the 
cubic nonlinear terms with n = 1 can not be treated by the nonlinear analysis in Strichartz 
spaces. Additional constraints must be imposed to ensure that small initial data in H 1 
scatter to zero. 

A similar question has been addressed in the context of the nonlinear Klein-Gordon 
equation, 

u tt -u xx + u + \u\ p ~ 1 u = 0, (5.1) 
Let us consider the semi-group of the linear Klein-Gordon equation, 

S(t) := e - lt{d *\ t>0, 

where (a?) = Vl + x 2 . The L°° - L 1 norm of the semi-group decays like O^ 1 / 2 ) as t — > oo. 
As a result, the term = 0(t _ ^ p_1 ^ 2 ) is absolutely integrable in t for p > 3. The 

dispersive decay of small initial data for p > 3 was proven by Georgiev & Lecente [5]. In 
the critical case p = 3 of the cubic nonlinear terms, the decay of small initial data and the 
scattering to zero was recently obtained by Hayashi & Naumkin [9] . 
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We will show that the results of Hayashi & Naumkin [8j [9] can be equally applied to 
the nonlinear Dirac equations (|2.ip . Using the Fourier transform, we rewrite the system as 



iu t — ku + v = f, 
ivt + kv + u = g, 



(5.2) 



Using the projection matrix 
P = 



1 -Vl + k 2 - k 

Vl + k 2 + k 1 



P- 1 



we can write (|5.2p in the equivalent form 



Vl + k 2 - k 1 

-1 y/lTW -k 



ia t - (k)a = ^{-f + ((k) - k)g), 
ibt + (k)b = j^(g + ((k)-k)f), 



(5.3) 



In the physical space, this system takes the form 



at + i(id x )a = %>(id x ) l (-f(a,b) + ((id x ) + id x )g(a,b)), 
bt - i{id x )b = ^d x )-\g{a,b) + ((idx) + id x )f(a,b)), 



(5.4) 



where 



f(a, b) = f(b+{id x - (id x ))a,a + ({id*) - id x )b), 
g{a,b) = g(b + (id x -(id x ))a,a + ((id x ) -id x )b). 

If / and g are homogeneous cubic polynomials in variables u and v, then / and g are 
cubic polynomials in variables a, b, d x a, d x b, {id x )a, and (id x )b. By Remark 1.1 in [9], all 
these cubic coefficients can be treated in spaces i?^(R) equipped with the norm, 

Hffr := \\(x) s (id x ) m u\\ L 2. 

The method of Hayashi &; Naumkin [8, 9j gives the following theorem. 

Theorem 5.1 Fix small e > and assume that uo G Hf(M) with ||uo||^4 < e. There exists 
eo > such that for all e £ (0, eo), the Cauchy problem for the nonlinear Dirac equations 
\2. ty) with the quartic function W in \2. S\) has a unique global solution 



u(i) 6 C(R+,Hf[ 
satisfying the time decay estimate 

\\(id x )u(t)\\ L o. <Ce(l + t)-^ 2 , ie 
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Remark 5.2 Small Hf norm on the initial data implies small H 1 and L 1 norms of\iQ. 

We conclude the article by considering the integrable case of the nonlinear Dirac equa- 
tions (|2.ip with W = 2\u\ 2 \v\ 2 , which is referred to as the massive Thirring model (MTM). 
Theorem 13,11 implies global existence of solutions of (MTM) in H . Theorem 15,11 implies 
that these solution scatter to zero in the L°° norm. More results were obtained for the 
massive Thirring model recently. 

Selberg and Tesfahun [2D] proved local well-posedness of (MTM) in H S (M.) for s > 
and global well-posedness in H S (M.) for s > \. Machihara et al. [14] proved for a similar 
nonlinear Dirac equations with quadratic nonlinear terms that local well-posedness holds 
in H S (M) for s > — \ and that the Cauchy problem is ill-posed in iif _1//2 (IR). Using ideas 
from [20] and |14j . Candy [2] proved local and global well-posedness of (MTM) in I? . 

In characteristic coordinates, 

x — t x + t 



2 1 2 1 

the massive Thirring model with W = 2|it| 2 |v| 2 is written explicitly by 

iu T + v = 2\v\ 2 u, 
—iv£ + u = 2\u\ 2 v. 



(5.5) 



Let us introduce the change of variables, 

u(£,t) = ~™(£,t) exp^- l -j^\w\ 2 (e,r)d^y (5.6) 
The second equation of system ()5.5j) can be solved with 

v&t) = -^^(^^exp^-^y 00 ^! 2 ^,^'), (5.7) 



where 



/•oo 



If v(-,t) 6 -fT 1 (R), then the zero-mass constraint J K u;(^,r)ci^ = has to be added. 

With the substitutions (|5.6p - (15.7p to (|5.5p . the massive Thirring model becomes the 
scalar evolution equation 

w T - d^w + ^l^ 1 ^! 2 ^ = 0. (5.8) 
The scalar equation (|5.8f) is invariant under the following change of variables, 

w = 5W(X,T), X = 5 2 £, T = 5~ 2 t, 5 > 0, (5.9) 

which implies that the massive Thirring model is the L 2 -critical model with ||w||x,2 = ||I / F||£2. 
Therefore, it is natural to expect that the dispersive decay to zero can occur already for a 
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smooth initial data with a small L 2 -norm. To deal with this question, we shall review the 
inverse scattering transform method for the massive Thirring model. 

The scalar equation in characteristic coordinates (|5.8p appears as a solvability condition 
[TH [12] of the spectral problem 



and the linear time-evolution problem 



" fa ' 




—iX 2 Xw 




' fa ' 


_ 1p2 _ 




—Xw iX 2 




. fa _ 



Or 



fa 


= i 


. fa _ 





W w \ 2 

rjd^ 1 w 



w 



w 



fa 
V'2 



(5.10) 



(5.11) 



where A G C is a (£, r)-independent spectral parameter and rj = Note that the mas- 
sive Thirring model in the laboratory coordinates can also be represented as a solvability 
condition of the Lax system (TUJ E3] . 

The spectral problem (|5.10p has some symmetries. If (fa(x; X),fa(x; A)) is a solution of 
(l5TT0j) . then 



(fa(x;-X),-fa(x;-X)), (fa(x;-X),fa(x;-X)), (fa (x; A), -fa (x; A)) 



(5.12) 



are also solutions of (15. 10H . 

The continuous spectrum of the spectral problem (|5.10p is located for A G R U {2R}, 
whereas isolated eigenvalues are located symmetrically in quartets (A, —A, A, —A) in all quad- 
rants of the complex plane for A [12j. Isolated eigenvalues are associated with solitons that 
occur in the long-time dynamics of the solution w(£, r) thanks to the independence of A 
from t and the inverse scattering transform technique. We shall prove that solitons are 
absent if w(-,t) for a frozen r has a small norm in L 2 (R) nL°°(R) (or if^R)). For clarity 
of presentation, we do not write r in the arguments of w and 01,2 • 

Lemma 5.3 Fix small e > and assume that w G L 2 (R) with \\w\\i2 < e. There is C > 
such that the spectral problem \5. 10\) admits no solutions in L 2 (R) for any A G C with 
arg(A) G (Ce 2 , § - Ce 2 ). If in addition, w G L l {R) n L°°(R) and d^w G L^R) with 

IMIzaGMU 00 + II^HIlO < e > (5.13) 

then the spectral problem \5. 1 0\) admits no solutions in L 2 (R) for any A G C. 

Proof. Let us choose A G C in the first quadrant of the complex plane. Because of the 
symmetry (|5.12p of eigenvectors, the results are valid in all four quadrants of the complex 
plane. If A is in the first quadrant of C, then Im(A 2 ) > and we can parameterize A by 
A = |A|e ie with 6 G (0,f). 

If Im(A 2 ) > 0, then e~ lX2 ^ decays to zero as ^ — > — oo, and we can introduce "*/>(£) = 
e~ lX ^if(^) with boundary conditions lim^-oo <£>(£) = (1, 0) T for eigenvectors of the spectral 
problem (|5.10p . 



11 



Integrating the system from — oo to £ under the boundary conditions for <p(0> we obtain 
the integral equations, 

MO = 1 + A [ wiOMC'H', MO = -A f e 2lX2 ^w(0)M0) d 0- (5.14) 

J —oo J —oo 

Using the exact integral 

/(A):=|A 2 | I" e- 2 ^>dx |Ah ' ' 



2Im(A 2 ) 2sin(20)' 

the Schwarz inequality, and Young's inequality for convolution integrals, we obtain 

|| - 1||l°° < II^IUallA^H^, ||A99 2 ||l2 < I(X) \\w\\ L 2 ||</?i||l°°- 
Closing the inequalities and using fixed-point arguments, we can see that if 

/(A)Hl£ 2 <i, 

then there is a unique solution of system (I5.14p for tp\ G L°°(R) and A922 G L 2 {E) such that 
II Pi — l||i°° < 1- Therefore, 921 (£) ^> as £ — >• +00, and so ~4>i(0 grows exponentially as 
£ — > +00. Eigenvectors in I? may only exist if J(A)||u> ||| 2 > U that is, if either 9 G (0, Ce 2 ) 
or 6 G (f - Ce 2 , 1) for some C > 0. 

To eliminate eigenvectors everywhere in the first quadrant of C, we add now the condi- 
tion (|5.13p . Integrating the second equation of system (|5.14[) by parts and using the first 
equation, we obtain 

XMO = "A 2 f e 2 ^-^w{0)M0)d0 
J —00 

1^)^(0-1 / ? e 2lX2 ^(d^w{i')M0)+w{0)d il M0))d0 



OO 



Using Holder's inequality and Young's inequality for convolution integrals, we obtain 
||yi - l||i°° < H^llLillA^IlL 00 , 

I|a^2||l°° < 2 (Iklli^ + II^^IIlO ll^ilk°° + ^IkllL 2 !!^!!^- 

If ||^|| 2 2 < 2, then there is C > such that 

||^i - 1||l°° < C||w|| L i(||t(;||L°° + H^-wlliOH^illLoo, 

II Ac^ 2 |U°° < c(l|w|U°° + II^^IIlOII^iIU 00 - 

Under the condition (j5. 13|) . there is a unique solution of system (|5.14p for cp% G L°°(M) and 
Xf2 G L°°(M) such that \\ip\ — < 1- Repeating the arguments above, we conclude the 

proof that no eigenvector in 1? exists for any A G C under the condition (|5.13p . □ 
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Remark 5.4 Using the scaling transformation A5.9\) . we can see that both \\w\\l2 and 
\\w\\li (\\w\\l°° + ||^w||^i) are invariant with respect to parameter 5. Soliton solutions of 
the massive Thirring model are supported by particular values for these quantities. IfWwW^ 
and \\w\\ L i (||u>||l°° + Hc^tuH^i) are below these particular values, no solitons can occur in 
the long-time evolution of the massive Thirring model. 

Further analysis of the inverse scattering transform using the time-evolution problem 
(|5.11j) may give an analogue of Theorem 15. II for the massive Thirring model, perhaps, with 
relaxed assumptions on the initial data vlq. Another interesting open problem is to explore 
global existence of the massive Thirring model in L 2 [2] and obtain L 2 -orbital stability of 
MTM solitons. A similar task was recently achieved by Mizumachi & Pelinovsky [TS] in 
the context of the nonlinear Schrodinger equation. The massive Thirring model is more 
interesting for orbital stability analysis of solitary waves. Because it is associated with 
the sign- indefinite Hamiltonian function (|3.ip . no orbital stability in H 1 can be extracted 
from the standard energy analysis. These open problems will likely to attract interests of 
researchers in near future. 
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